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Extended BTK-type theoretical methods were used for the calculation of 
conductance spectra of NS and NIS junctions, where S is either a d-wave or 
a d+is-wave superconductor, and N is a normal metal. We found that the 
length scale over which the bulk d-wave order parameter varies spatially near the 
interface (the "decay length") has a significant effect on the main features of 
the conductance spectra. A relatively large decay length of the order of 10-20£o 
together with a small zs-wave component, can explain many peculiarities observed 
in the experimental conductance spectra of various junctions. We attribute the 
long decay length to either a reversed proximity effect in the high transparency 
junctions, or to the presence of a low T c < 10K s-wave component with a much 
longer £ 0s ~ 10£od near the interface. 
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INTRODUCTION 

Many investigations of the special properties of d-wave superconductors and 
their influence on the conductance spectra have been reported in recent years. |U 
2J A large number of them deal with NIS or SIS junctions in the tunneling limit 
where the effective BTK barrier strength 

is larger than ~1 (strong barrier). [3J Here H is the amplitude of the 5-function 
barrier at the interface. From a theoretical point of view, the tunneling limit is 
simpler to treat because in this limit of weak coupling between the electrodes, 
one can consider the electrodes as almost independent entities. More specifically, 
one can assume that the spatial changes of the superconducting order parame- 
ter, arising from the presence of a surface, occur only in a very small distance 
from the interface inside the superconductor. Thus, to a first order approxi- 
mation these changes are often ignored. Conductance spectra of NIS junctions 
where S is a d-wave superconductor, are characterized by a zero bias conduc- 
tance peak (ZBCP) in the node-direction, and by a gap in the anti-node (main 
axis) direction. Theoretical calculations of low transparency junctions reproduce 
these features quite well, also in the simpler case when a spatial independent 
pair potential is taken into account. |4J The ZBCP is a manifestation of zero- 
energy bound states of quasi-particles near the interface. The appearance of 
these bound states along the node direction of a d-wave superconductor results 
from the sign change of the energy gap A(<p + ) = — A(^?_), where tp± are the 
angles between the wave-vector of the incoming or scattered quasi-particle and 
the normal to the interface, and cp + = n — ip^. 



Theoretical calculations of the conductance spectra of junctions with a weak- 
barrier are more problematic. The problem arises from the strong mutual depen- 
dence of the pair potential on both sides of the NS interface. This is in contrast 
to the strong barrier case, where the two electrodes are almost decoupled. Thus, 
in order to obtain an explicit form of the spatial dependence of the pair potential 
(SDPP), a complete knowledge of the quasi-classical Green function is necessary, 
not only in the S side but also in the N side of the junction. Therefore, a proximity 
effect must be directly included in the calculations. There is no consensus on the 
peculiarities of the proximity effect in junctions with anisotropic superconductors 
at the present time. The penetration length of the order parameter into the 
N side and its symmetry there, depend sensitively on the specific properties of 
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the normal metal. [5, HO |H1 IS] Recently, it was found, that even ferromagnetic 
or antiferromagnetic metals in NS junctions have a long-range proximity effect 
under certain conditions, [10] HU H2] and this is in contrast to the predictions 
of certain theories. [13] For the 123-family of superconductors, the length scale 
of the proximity effect in the N side depends significantly on the oxygen-doping 
level of the S side. This results from the fact that in underdoped YBCO T c is 
reduced compared to optimal doping, and the value of the coherence length is 
increased. A-priori, the proximity interaction between the two electrodes should 
depend strongly on the barrier strength. But as was shown previously. [141 IT51 Ho] 
the effect of the barrier is not always so evident, and even a barrier of finite 
strength is not always effective in reducing the conductance of the junctions. In 
particular, in NS junctions with similar materials and similar density of states 
(DOS) on both sides, the influence of a barrier is suppressed significantly. This 
occurs when both electrodes are made of differently doped cuprate compounds. 
According to Deutscher and De Gennes,[17j in such cases the amplitude of the 
order parameter in the normal side near the interface can reach a value close 
to that in the superconductor, and there is almost no discontinuity or jump at 
the interface. Therefore, in this type of junctions the diffusion of pairs from the 
superconductor into the metal is almost unaffected by the barrier, and the decay 
distance of pairs in the normal phase may be anomalously large. 

YBCO has a short coherence length £ = Jivf/Aq of about 10 — 30A where 
A is the bulk value of the d-wave pair potential . [18J In spite of this fact, it follows 
from the above discussion that in studies of NS-junctions with similar cuprates 
one should examine a wide range of decay lengths in S and proximity lengths 
in N. To clarify this different length scales dependence, we decided to calculate 
the corresponding conductance spectra in a phenomenological way. We didn't 
develop a new model for self-consistent calculations of the pair potential in prox- 
imity systems. Instead, we assumed reasonable pair potential shapes (SDPP), 
and calculated the resulting conductance spectra. From the results we were able 
to identify which of the peculiarities of the SDPP had the most significant effect 
on the conductance spectra, and what form of SDPP leads to the best fit of the 
experimental data. 

RESULT AND DISCUSSION 

We have calculated numerically the conductance spectra of a normal metal 
(N) - d-wave superconductor (S) junction as a function of applied voltage us- 
ing the basic model developed by Tanaka and Kashiwaya.[lJ According to their 
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theory, the conductance is given by: 



tt/2 oo 

/ / cr s (E, v )(-?li§t?4-)a N (v)co S <pdvdE 

««) = ^ ~ 2 , (2) 

f <7jv(<£>) cos tpdip 
-tt/2 

where the integrations are over the angles between the trajectory of the incoming 
quasi-particle and the normal to the interface of the junction, and the energy. 
The normal metal conductance a^, characterizes the transparency of the junction 
and is given by: 

4cosV /q\ 

aN = ~A 2 T^2- ( S > 

4 cos^ if + Z 2 

where Z is the effective barrier strength as given by Eq. (1). f(E) in Eq. (2) 
is the Fermi distribution function at temperature T. The function as in Eq. (2) 
represents an extension of the standard BTK expression for the conductance |3J 

a s = l-\b\ 2 +\a\ 2 , (4) 

to the case of spatial dependent order parameter. In Eq. (4), b and a are the 
coefficients of the normal and Andreev reflections, respectively. In the simple 
BTK model, the reflection coefficients are expressed in terms of the functions 
r±(E,(p) in the following way: 

^5 = l + a N \T + \ 2 + (cr N -l)\T + T_\ 2 

<? N |i + K-i)r + r_| 2 ' U 

The r± in this equation are functions of the energy E and angle ip only, since the 
pair potential is assumed to be spatially constant. In the present study however, 
we use the extended BTK model which takes into account the spatial dependence 
of the superconducting pair potential. Following Kashiwaya and Tanaka,[lJ we 
calculated the F± values at the interface by solving a set of Riccati type equa- 
tions (Eqs. (3.50 - 3.52) in Ref. [lj) in which the spatial dependent pair potentials 
A±(ip,x) appear as coefficients (x is the distance from the interface in S). 

In our simulations we have chosen the following functions to represent the 
spatial dependence of the pair potential of the dominant d x 2_ y 2 and sub-dominant 
s components of the order parameter: 
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: = tanh((C + A; 1 )*A;2), (6) 



and 

X( ° -A; 3 exp(-A; 4 *C 2 ), (7) 



A 

where ( = x/£ is a dimensionless length (see Fig. 1) . The superconductor oc- 
cupies x > and the interface is exactly at x=0. Here the ki (i = 1..4) are 
dimensionless parameters which depend on the junction orientation relative to 
the crystalline axis. The modelling functions in Eqs. (6) and (7), satisfy the main 
boundary conditions on both component of the pair potential (A<^ 2 2 — > A 
and A s — > as x — > oo), and are consistent with the self consistent model 
results. [2UJ [Tj5] Moreover, we note that the fine details of the shape of the SDPP 
is not essential, as it does not affect the resulting spectra in any major way. 
In order to confirm this statement, we examined different possible functions by 
replacing, for instance, tanh(£) with 1 — exp(— Q. Results of these calculations 
show that such changes lead only to minor changes in the energies of the high- 
order bound states in the conductance spectra. We found out that the most 
important parameters of the SDPP shape which affect the conductance spectra, 
are the maximal amplitude of the sub-dominant component (the k 3 coefficient 
in Eq. (7)) and the quasi-particle decay length (QPDL), which represents the 
characteristic length scale over which the pair potential varies spatially near the 
interface (governed by ki,k 2 and k± in Eqs. (6) and (7)). 

In all our calculations we have taken into account the thermal smearing factor 
which is due to the finite (low) temperature, and the finite life time broadening 
of the quasi-particles at the Fermi surface. It is a well known fact from the theory 
of Fermi liquids, [21J that the inverse lifetime of the quasi-particles is proportional 
to E 2 . In real superconductors, there are some quasi-particle relaxation processes 
which include electron-electron and electron-phonon interactions, which modify 
this simple square low dependence. Kaplan et a/.|22j showed that the general 
dependence of such relaxation rates can have a very complicated dependence on 
energy which is unique for each material. Dynes et al. however, used successfully 
a constant lifetime broadening in the theoretical fits of their measured conduc- 
tance spectra in low T c SIS junctions. |23J Such approximation indeed works well 
for fitting of experimental spectra in isotropic superconductors and in the low 
transparency limit. In this case, A(x) ~ constant, the conductance is propor- 
tional to the quasi-particle density of states, and this has strong sharp peaks at 
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E = ±A . The addition of lifetime broadening as a small imaginary part to 
the energy allows one to avoid divergence at the peak energies and does not 
produce any essential change in other regions. Our attempt however, to apply 
this method to high-transparency spectra was unsuccessful. We attribute this 
to the appearance of additional bound states due to a significant increase in the 
decay length (see below). It is thus clear that in such cases, taking into account 
the energy dependence of the lifetime broadening is essential. In our calculations 
we used a fourth degree polynomial approximation for this dependence: 



where the coefficients 7, < 1 (in units of A ) were chosen in order to obtain the 
best fit to the experimental conductance spectra. Generally, the function T is 
added as an imaginary part to the energy in the expression for the DOS. [23] H] 
In our calculations however, the equivalent of the DOS is the function a s of 
Eq. (2). Since we obtain as by solving two coupled differential equations in 
which the energy is a parameter, we could not just add the iT term to E directly. 
The fact that the boundary conditions of these equations include nonlinear terms 
in E, complicates things further. We therefore calculated as first without any 
broadening term, and then added the broadening numerically by a convolution 
of the obtained as(E) function for each angle (p, with a Lorentzian distribution 
function: 



of width T(E). 

We focus now on the case of a pure d-wave symmetry of the order parameter 
without any sub-dominant component. Two basically different cases of the node 
(a = vr/4) and anti-node (a = 0) orientations will be discussed, a is the angle 
between the lobe-direction of the d x 2_ y 2-wave order parameter and the normal 
to the interface. In the node direction, a ZBCP exists independently of the 
values of the coefficients k\ and k 2 in Eq. (6), and the effective barrier strength 
Z. Nevertheless, by varying k\ and k 2 in Eq. (6), the shape of the calculated 
conductance spectra changes. Decreasing their values, leads to an increased 
suppression of the order parameter at the interface and an increased decay length. 
As a consequence, additional quasi-particle bound-states are created. [24 1 [T^J 
These new bound-states are formed in a potential well defined by the interface 
on one side and the growing d-wave pair potential toward the bulk value on the 



T{E) = 7l E 2 + 72 E 3 + 



(8) 



L(E) 



1 T(E) 
n(E-E o y + T\E) 



(9) 
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other side. They manifest themselves as additional peaks in the function o- s {E) 
versus E. Their number and energies vary for different angles if, and according 
to Matsumoto and Shiba,[19j their energies are given by 



Tl Tl 

E n = ±A | sin(2</?) \ \ -(2 ), (10) 

1 ' ty rp 



where n is an integer, 

2 | sin(2y>) 



(11) 

and < n < r. k 2 in Eq. (11) is different from zero since a k 2 =0 in Eq. (6) 
would yield A = which means that no superconductivity exists. The number 
of bound states is determined by r. The larger the value of r the larger the 
number of these states. Eq. (11) does not give the dependence of r on fcj, 
because in Ref. [19J only the special case of a superconductor-vacuum interface 
is considered (Z — ► oo). Therefore in Ref. [19J, k\ is always zero in the node 
direction. It is however clear from Eq. (6) that with increasing k\ the potential 
well is filled up, the spatial dependence of A approaches that of a step function, 
and the bound-states disappear. Thus, in the general case of arbitrary Z and a, 
the number of bound-states is proportional to the depth of the potential well. 
This depends on both coefficients k\ and k 2 , and yields a more general expression 
for r as follows: 

r = 2|s y 2y?)l [l-tanh(fc 1 fc 2 )]. (12) 
k 2 

Eq. (12) shows that while we consider quasi-particles within a narrow cone 
around <f = (like in SIN junctions), bound-states will play a less significant 
role in the conductance (if the cone angle if — > then there are no bound 
states since r — > 0). In the weak barrier limit though, a broad cone of about 
— 7r/2 < ip < 7r/2 should be used, and the influence of bound-states on the 
conductance spectra becomes substantial as seen in Fig. 2. On each curve of 
this figure only two symmetrical maxima besides the ZBCP are present. The 
curves in Fig. 2 however, show the a(eV) dependence on V and not the &s{E) 
dependence on E. In the later, for each value of f a different number of addi- 
tional maxima exist. After integration of as over all angles ip, only two main 
peaks versus energy are left which contribute to the resulting maxima Epd in 
a(eV) as seen in Fig. 2. These side band peaks are the second most prominent 
features in the spectra. With increasing r, the voltage of these peaks shifts to 
low bias. Because of this, after the integration in Eq. (2) the ZBCP becomes 
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somewhat smaller and its width decreases. It is important to note, that in spite of 
the shift of additional maxima towards zero bias with increasing decay length of 
the pair potential, they can't annihilate the ZBCP, since the zero-energy bound 
state is robust and exists for all values of cp, while the other bound-states are 
not. Therefore, by changing the parameters k\ and k 2 we can't obtain splitting 
of the ZBCP in the node direction in the pure d-wave case. 

Splitting of the ZBCP along the node direction was however measured in SNS 
junctions. |25j In order to fit this data, it was necessary to add a sub-dominant 
component to order parameter with a different symmetry. Because of previous 
results, [22J I2Z] we nave choosen for this purpose to add an imaginary s-wave 
component to the pair potential. Now the form of the resulting order parameter 
depends on all four coefficient ki,k 2 ,k 3 and fc 4 . Our calculation results along 
the node direction and with a small Z value are shown in Fig. 3. The main con- 
clusion from this figure is that the sub-dominant component can not essentially 
influence the conductance spectra as long as the decay length of the dominant 
d-wave component is ~ £ - This behavior is different from that observed in the 
strong barrier limit, where a small decay length for both components is enough 
for obtaining splitting of the ZBCP (not shown). Tanuma et al. have obtained 
no sub-dominant component for low values of Z (Z ~ 1),[20J even though their 
Z value is defined as twice the BTK Z value that we use here. In addition, the 
order parameter that they obtained, almost always had a decay length of order 
~ £q. However, when we repeated their self consistent calculations for high 
transparency junctions, we could find cases (when T cs = 0.5T c d for instance) in 
which the s-component did exist and had a long penetration length of the order 
of 10£ into the bulk. Moreover, in their self-consistent calculations of the order 
parameter they did not take into account the proximity effect. Therefore, the 
absence of a sub-dominant component in the results of Tanuma et al. for weak 
barriers, is not contradictory to our results. As can be seen in Fig. 3, we obtain 
the splitting effect by the s-wave component only when d-wave decay lengths 
longer than a few £ are used. 

A sub-dominant is component splits the ZBCP into two peaks, whose energy 
\E ps \ depends on the values of k 3 and k^ in Eq. (7). There is no exact analytic 
formula for this dependence, but numerical calculations yield an approximate 
relation: 

\E ps \ ~ Cifc 3 exp(-c 2 /c4) for k 3 < A , (13) 
where c\ and c 2 are constants which depend on Z. When two components of the 
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order parameter coexist, the conductance spectra is a superposition of contribu- 
tions containing the bound-states of both. These bound-states can annihilate 
one another unless the value of E ps is close to the energy E p d of the bound- 
states of the dominant component. It thus follows that the smaller the s-wave 
component, the larger the required decay length of the d-component, which is 
needed in order to obtain a clear splitting. 

A simulation of a the conductance spectra fitting our previous experimental 
data along the node direction is given in Fig. 4. The original data was obtained 
for SNS junctions. [25J Since no theoretical model is available at the present time 
for the calculation of the conductance of SNS junctions which takes into account 
the finite decay length of the pair potential at the interface, we used as a first 
approximation the SN model discussed here, but with a re-normalized energy 
scale to fit the SNS data. This re-normalization was simply done by doubling 
the energy scale in the SN calculations (multiplying it by a factor of 2). We can 
thus compare the modified SN simulation results with the measured SNS data 
on the same energy scale of E/Aq. Coexistence of both d and is components 
of the order parameter near the interface had to be assumed in order to fit the 
data. The good agreement between theory and experiment indicates that large 
decay lengths of the d-wave pair potential are actually present near the interface 
of real junctions. 

We now turn to junctions along the anti-node direction. This direction is 
characterized in the tunneling regime by a "V" shaped conductance spectra with 
a minimum at zero bias whose depth is a function of the barrier strength Z. 
For a low transparency barrier a(E = 0) ~ 0, but as Z decreases this value 
increases. For weak barriers (small Z) the conductance is controlled mostly by 
Andreev reflections, and a(E = 0) reaches a value of 2 when 6 = and a = 1 in 
Eq. (4). Again, we varied the values of the coefficients fcj (i = 1..4) in order to 
study the influence of the order parameter near the interface on the conductance 
spectra. Qualitatively, all our previous reasoning remain valid. Decreasing ki 
and k 2 leads to the appearance of bound-states which are seen as new maxima 
in the simulated conductance curves. Similar to the results along the node direc- 
tion, we find that with increasing decay length of the d-wave pair potential, the 
bound-states appear at lower bias. But they never occur at zero-bias, where the 
conductance minimum persists up to very low Z values. Fig. 5 shows simulations 
of conductance spectra for a weak barrier of Z = 0.5 along the (100) direction 
for different values of the coefficient k 2 while k\ = 0. The main feature of these 
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spectra is that the peaks at the sub-gap energies increase with decreasing k 2 
(increasing decay length). Moreover, addition of a small sub-dominant s-wave 
component to the d-wave order parameter when a large decay length is used, 
does not practically change the form of the conductance spectra (The solid and 
dotted curves in Fig. 5 are almost indistinguishable). This leads to an important 
conclusion that the conductance along the (100) direction in high transparency 
junction with a long-range proximity effect, is basically insensitive to the presence 
of possible sub-dominant components of order parameter. 

Fig. 6 shows experimental conductance spectra of a low-transparency junc- 
tion along the (100) direction,[25j together with a best fit simulation for this 
case. The energy scale of the SIN simulation was re-normalized to the equiva- 
lent SIS scale of the experiment as before (was multiplied by a factor of 2). The 
low transparency is a result of a large normal resistance R N which is found in 
the experiment for the (100) orientation. The roughness of the interface in this 
type of junction was found to be smaller than the coherence length £ . thus no 
faceting can lead to the low transparency. Twinning however, can cause differ- 
ent Rn values in different directions. We found out that a pure d-wave order 
parameter neither with a small nor with large decay length can fit our data well, 
especially not at low bias (see the dash and dotted curves in Fig. 6). It turned 
out that the contribution of bound states is essential in order to describe cor- 
rectly the low energy dependence of the experimental conductance spectra. A 
typical contribution of bound states is shown in Fig. 6 by the dotted curve for 
a 20£o decay length. For the best fit in this figure we used a superposition of 
mainly a constant d-wave component (76%), and the additional contributions 
from three longer decay lengths of 3£ , 10£o ar| d 20£ (8% of each), and found 
a remarkable agreement between the simulation and the experimental data (the 
solid line in Fig. 6). The fact that we had to use such a superposition means 
that in reality, at least in the lobe direction, there is some non-uniformity in the 
values of the decay length over the junction cross section. This can be due to 
the twinning mentioned above, to a non-uniform oxygen distribution or to some 
other disturbances. We note that in grain boundary junctions a nonuniform dis- 
tribution of different barrier width was found, [28J which in a way is similar to the 
decay length distribution observed here. 

In the following we summarize the main properties of SN and SIN junctions 
resulting from our simulations: 

i)For a pure d-wave superconductor in the node-direction, large decay lengths 
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of the order parameter near the interface lead to additional peaks in the con- 
ductance spectra at sub-gap biases. The peak energies shift to lower bias with 
increasing decay length of the order parameter, while the ZBCP amplitude and 
its width decrease. 

ii) ln order to obtain a well defined splitting of ZBCP in the conductance spectra 
along the node direction as seen in the experiment, one has to add a small sub- 
dominant is component to the order parameter, and keep a long decay length. 

iii) ln contrast, using a pure d-wave with a large decay length in high transparency 
junctions along the lobe direction, the gap peaks in the conductance move to 
lower energies and this looks like splitting of the broad Andreev-enhanced con- 
ductance regime at E < A . 

iv) The energy dependence of the quasi-particles lifetime broadening has to be 
taken into account in the simulations of the conductance spectra. 

v) The fact that a long decay length in the superconductor had to be used in 
order to obtain a good agreement between the theoretical simulations and the 
measured conductance spectra, gives further support to the notion that a long- 
range proximity effect exists in the cuprates. 

CONCLUSIONS 

In the present study we investigated how the spatial dependence of the order 
parameter near the interface of junctions affects the resulting conductance spec- 
tra. We found that the decay length of the d-wave pair potential is the most 
significant parameter affecting the spectra. More specifically, for high trans- 
parency junctions along the node direction, the use of a long decay length of 
the order of 10 — 20£ yields a splitting of the ZBCP, provided an additional 
small is component exists near the interface. If in contrast a ~ £ decay length 
is used, no such splitting develops in the simulations, and this is in contradic- 
tion to experimental observations. For low transparency junction along the main 
lobe direction, a small decay length of the order of ~ £ could not describe the 
data neither with a pure d-wave nor with a d+is order parameter. Only when a 
distribution of decay lengths between 0£ and 20£ was assumed, contributions 
from additional bound states led to a good fit of the experimental data. Another 
important parameter that affects the conductance spectra strongly is the quasi- 
particles lifetime broadening and its non-trivial energy dependence (F(E)). We 
actually modified its basic Fermi liquid E 2 dependence to fit the experimental 
data along the node direction, and then used it successfully without any further 
changes to fit the lobe data. Finally, we note that a long decay length in the S 
side is consistent with a reversed proximity effect in high transparency SN junc- 
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tions, and the low T c of the is component. 
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Figure Captions 



Fig. 1: Typical spatial dependencies of the d-wave and s-wave pair potentials 
near the interface, used in our simulations. The pair potential axis is 
normalize to A , the bulk value of the superconductor's gap. 

Fig. 2: (color) Simulations of normalized conductance spectra of a high trans- 
parency SN junction with Z=0.5 along the node direction for a pure d-wave 
symmetry of the pair potential. Different curves correspond to different 
values of the the decay length. 

Fig. 3: (color). Theoretical curves of normalized conductance spectra of a high- 
transparency NS junction along the node direction using a d+is symmetry 
of the pair potential. In each curve the decay length of both components 
of the pair potential is the same, except for the dash double dotted line 
where the decay length of the d-wave component is £ while that of the 
s-wave component is 18£ - 

Fig. 4: (color) Measured normalized conductance spectra of a high transparency 
YBa2Cu 3 O e+x /YBa2Fe 0A5 Cu2.55O e+x /YBa 2 Cu 3 O 6+x junction along the 
node direction (symbols) together with a best fit simulation (line). The 
broadening T versus E is given in the inset. 

Fig. 5: (color) Calculated normalized conductance spectra of a high-transparency 
SN junctions with Z=0.5 along the lobe direction. All curves were obtained 
by the use of a pure d-wave symmetry of the pair potential with various 
decay length, except for the dotted line where a d+is symmetry was used. 

Fig. 6: (color) Measured normalized conductance spectra along the lobe direc- 
tion (symbols) of a low transparency junction of 

YBa 2 Cu 3 OG +x /YBa2Fe 0A5 Cu2.r J 5OQ +x /YBa 2 Cu 3 O 6+x together with a 
best theoretical fit (solid line). The energy dependent lifetime broadening 
T(E) is the same as that in Fig. 4. Also shown are the conductance con- 
tributions of a constant d-wave (dash line) and a d-wave with 20£ decay 
length (dotted line). 
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